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Abstract. Three-gluon to three-gluon scatterings lead to rapid thermahzation 
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1. Introduction 

While the energy of nucleus-nucleus collisions increases, the density of initially cre- 
ated partons gets larger and larger. Number density of gluon matter has been 
required to be more than 30 fm~^ in RHIC Au-Au collisions and may reach 140 
fm~"^ in LHC Pb-Pb collisions. At such high gluon number densities, new physics 
processes begin to play an important role and are responsible for some phenomena. 
Three-gluon to three-gluon scatterings have been shown to lead to rapid thermaliza- 
tion of gluon matter created in RHIC central Au-Au collisions [[T]. Thermalization 
has been studied from parton cascade models El U and in other model attempts 

[El El El m n cni ini m 

To study triple-gluon elastic scatterings, we need to know how frequently the 
triple-gluon scatterings occur. This is shown by the ratio of three-gluon to two-gluon 
scattering numbers in the next section. The triple-gluon and triple-quark elastic 
scatterings and transport equations are given briefly in Sections 3-5, respectively. 
Summary is in the last section. 
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2. Ratio of Scattering Numbers 

we estimate the ratio of numbers of scatterings occurring between two gluons or 
among three gluons at the moment when gluon matter is formed. The counting of 
scatterings is made for gluons in an anisotropic momentum distribution obtained 
from HIJING simulation [J^ for initial central Au-Au collisions at ^/snn = 200 



where Ra — 6.4 fm, tini ~ 0.2 fm/c and rapidity region | y |< 1^ = 5. 9{x) equals 
for a: < or 1 for a; > 0. One thousand gluons are generated from the distribution 
within the cylinder which has a radius of Ra and a longitudinal regime of -0.2 fm 
< z < 0.2 fm. Then the maxima of the numbers of the two-gluon and three-gluon 
scatterings are 500 and 333, respectively. We determine a scattering of two gluons 
when the distance of the two gluons is less than a given interaction range. If three 
gluons are within a sphere which center is the center-of-mass of the three gluons 
and which radius equals the given interaction range, a scattering of the three gluons 
occurs. It is plotted in Fig. 1 the ratio of the numbers of the three-gluon scatterings 
to the two-gluon scatterings at the time <ini. If the interaction range approaches 
zero from 0.1 fm, the number of the three-gluon scatterings reduces faster than the 
two-gluon scatterings; when the interaction range is larger than 0.6 fm, both of the 
scattering numbers are close to their maxima; when it is larger than 0.1 fm, the 
ratio varies around 0.7. The importance of the three-gluon scatterings is verified. 

3. Three-Gluon to Three-Gluon Scatterings 

The triple-gluon elastic scatterings involve many diagrams at order [ ^ and 
only two of them are selected to be shown in Fig. 2 for illustration. The scattering 
processes at tree level contain three-gluon and four-gluon couplings. Three incoming 
gluons interact at different space-time points or at the same space-time point. The 
squared four-momenta of propagators, (7^21 922 ^^'^ 923; tend to zero which 
causes Coulomb exchange divergence. The divergence is removed with the use of 
a screening mass. The triple-gluon scattering B_ cannot be thus identified as an 
iterative process of two successive scatterings of on-shell gluons. Such screening 
mass is evaluated in the use of a formula in Ref. [IT^. 

Squared amplitudes of three-gluon scattering diagrams are derived with For- 
tran code in the Feynman gauge. Interference terms of different diagrams are also 
calculated. If a gluon's four-momentum is labeled as pi = (£'i,Pi) in the pro- 
cess g(pi) + g{p2) + gips) g(p4) + g(P5) + giPe), squared amplitudes are out- 
puted in terms of the following Lorentz-invariant variables, S12 = {pi +P2)'^,S23 = 

{P2+P3f,S3i = {P3+Plf, Ui5 = (P1-P5)^,U16 = {Pl~P&f, "24 = {P2~PiY,U2& = 

{p2 — Pof' , U34 = (p3 — pnY and U35 = (pa — p^Y . In derivations, all possible 
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Fig. 1. Ratio of three-gluon scattering to two-gluon scattering numbers versus the 
interaction range. 

exchanges of final (initial) gluons are also taken into account to obtain different 
diagrams. 




B_ B 
Fig. 2. Scatterings of three gluons. 
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4. Three-Quark to Three-Quark Scatterings 

We use quark-quark elastic scatterings [ I15L and three-quark to three-quark 
elastic scatterings [El to study quark matter which is an ingredient of quark-gluon 
matter created in central Au-Au collisions at y'sjvw = 200 GeV. Quark matter is 
considered as consisting of only up-quarks and down-quarks and the two kinds of 
quarks have the same distribution functions. The triple-quark elastic scattering 
processes are plotted in Fig. 3. Exchanges of final (initial) quarks generate forty- 
two different diagrams . Squared amplitudes for all the diagrams are derived from 
Fortran code and expressed in terms of Si2, S23, S31, U15, uie, 1*24, U26, U34 and U35. 
The 42 diagrams contribute to the scatterings of three quarks with the same flavor. 
If one quark's flavor differs from the other two, only 14 diagrams give contributions. 



I I I I 



I I I I 




A_ Ajjt 
Fig. 3. Scatterings of three quarks. 



5. Transport Equations 



A transport equation of Boltzmann type including three-parton to three-parton 
scatterings is 



(Pp2 d^P3 d^Pi 



2-c/ig22 J (2TT)-^2h2 (27r)'^2£/3 (27rj'^2£/4 
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X I M2^2 P [/l/2(l ± /3)(1 ± /4) - fsMl ± /l)(l ± /2)] 

f d?P2 d^P3 <PPA d^P5 d^Pe 
2Eig33 J {2^Y2E2 {2Tr)^2E3 {2TTY2Ei {2^Y2E^ {2ttY2E^ 

X {2'Kf5'^{pi +P2+P3-P4~P5- Pe) I Ms^s |^ 

X [/i/2/3(l ± /4)(1 ± /5)(1 ± /e) - /4/5/6(l ± /l)(l ± /2)(1 ± h)] 

(2) 

where g is the color-spin degeneracy factor, the velocity for massless partons vi = 1, 
§22 = '^out' ^i^d g33 = ?^in!^out! wherc Hq^j (riout) is the number of identical final 
partons of 2 — > 2 (3 — *■ 3) scatterings and rtin for the 3 — > 3 scatterings is the 
number of identical initial partons except the parton in the distribution function 
/i- I A^2^2 P and | Mz^s P represent squared amplitudes of 2 ^ 2 and 3^3 
parton scatterings, respectively. The 3 — > 3 scattering processes involve a larger 
phase space than the 2^2 scattering processes. 
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Fig. 4. Gluon distribution functions versus momentum at tiso — 0.65 fm/c. The 
dotted, dashed and dot-dashed curves correspond to the angles 9 — 0°,45°,90°, 
respectively. The solid curve represents the thermal distribution function given by 
Eq. (3). 



Anisotropy of a parton momentum distribution as in Eq. 1^ can be eliminated 
by elastic scatterings among partons. The transport equation is solved from the 
time iini when anisotropic parton matter is formed and until the time tiso when 
local momentum isotropy is established. For gluon matter, gluon distribution func- 
tions at tiso — 0.65 fm/c in three different directions are shown in Fig. 4 by the 



6 



X.-M. Xu 



dotted, dashed and dot-dashed curves. The three directions correspond to the an- 
gles 9 — 0°, 45°, 90° with respect to an incoming gold beam direction. The solution 
of the transport equation at tiso exhibits similar momentum dependences in different 
directions relative to the beam direction. Such dependences can thus be fitted to 
the Jiittner distribution which differs substantially from the anisotropic momentum 
distribution given by Eq. 

where temperature T — 0.75 GeV, fugacity A — 0.065, a thermalization time of 
^iso — iini = 0.45 fm/c for gluon matter. When quark matter evolves independently, 
2^2 and 3^3 quark scatterings result in a thermal state with T = 0.59 GeV, 
A = 0.04 and tiso — ^ini ~ 1-8 fm/c for quark matter The thermalization times 

for gluon matter and quark matter are very different. 
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Fig. 5. Gluon distribution functions at tjni = 0.2 fm/c are shown by the dotted, 
dashed and dot-dashed curves which correspond to the angles 9 = 0°,45°,90°, 
respectively. The 2 — > 2 gluon scatterings result in the solid curve for the gluon 
distribution function at 9 = 0° at tiso = 0.65 fm/c. 

To definitely realize a role of the 3 — > 3 gluon scatterings, the gluon distribution 
function at the time ^iso = 0.65 fm/c resulted only from the 2 — > 2 gluon scatterings 
is shown as a solid curve in Fig. 5. The small discrepancy of the solid and the 
dotted curves exhibits a small variation of the gluon distribution function at the 
angle 6* = 0°. The gluon distribution functions at the two angles 9 = 45°, 90° at 
tiso — 0.65 fm/c are very close to the dashed and dot-dashed curves and are thus 
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not plotted. The variation of gluon distribution function resulted from the 2^2 
gluon scatterings is small in the forward direction and even negligible away from 
this direction. Therefore, in the thermalization driven by both the 2^2 and 3^3 
gluon scatterings, the gluon distribution variation from the anisotropic form to the 
thermal state is mainly caused by the 3^3 gluon scatterings. 

6. Conclusions 

We have studied thermalization of gluon matter and quark matter with the three- 
parton to three-parton elastic scatterings. The three- gluon elastic scattering pro- 
cesses give a considerably larger variation of the gluon distribution function than 
the two-gluon elastic scattering processes while the gluon number density is high. 
The triple-gluon scatterings are important at RHIC energies and yield the effect of 
rapid thermalization. The triple-quark scatterings give a variation of quark distribu- 
tion function comparable to the one that resulted from the quark-quark scatterings. 
Quark matter itself cannot thermalize rapidly at RHIC energies. 
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